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Abstract 

The effect of electron correlation on the photocurrent of self-assembled 
InAs/InGaAs quantum dot infrared photo-detector (QDIPs) is studied. It 
is found that Coulomb interaction and level mixing in the many-body open 
system lead to double peaks associated with the intra-band transitions involv- 
ing two lowest levels of the quantum dot. Furthermore, the photocurrent is a 
nonlinear function of the steady-state carrier density and it displays a plateau 
due to Coulomb blockade. 
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Recently, many efforts have been devoted to the understanding of transport properties in 
quantum dot photo-detectors (QDIPs)[l,2]. The advantage of the QDIP over quantum well 
photo-detectors (QWIPs) is that light can be directly coupled to the electrons in the normal 
incidence geometry due to the effect of QD confinement in directions perpendicular to the 
growth axis and the dark current is smaller for the same detection wave length considered [3]. 
Other significant features that are unique to QDs include the Coulomb blockade effect [2] and 
phonon bottleneck[4]. 

Due to the localized nature of electrons in QDs, it is essential to take into account the 
effects of Coulomb blockade in the analysis of photoresponse of QDIPs, which in general 
can be ignored in QWIPs. For the non-equilibrium system considered here, it is convenient 
to use the Keldysh Green function to calculate the transport and optical properties while 
including the electron correlation. This technique has been used extensively in the study of 
nonlinear transport properties of quantum systems [5, 6]. 

We solve the Anderson Hamiltonian of a two level system coupled with leads in the 
presence of an electromagnetic radiation with frequency uj. We adopt the approach given by 
Jouho et al[5,6] and extend it to the present case with asymmetric tunneling rates. We find 
that the time averaged tunneling current is given by 

(J(t)j) = (If - rj)^p _ e / ^[Tff L - Vff R \ImG^eM (1) 

where (Nj(t)) is the electron occupation number at QD, = f(e — /il) and fn = /(e — 
/jLr) are the Fermi distribution function of the left lead and right lead, respectively. The 
chemical potential difference between these two leads is related to the applied bias (V a ) via 
I^l — = eV a . Y 1 - and Yf denote the tunneling rates from the QD to the left and right 
leads, respectively, for electrons in level j. G r j(e,uj) is the retarded Green function for an 
electron in level j of the QD. The first term in Eq. (1) for the excited state (j = 2) provides 
the photo-induced tunneling current which exists only when 7^ , a condition that can 
occur in a system with asymmetric potential. 

Since the incident radiation considered in QDIP application is usually very weak, we can 
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ignore the renormalization of the retarded Green's function due to electron-photon interac- 
tion . Thus, we have (within the Hartree-Fock approximation) 



e-E 1 + iF 1 /2 e- E 1 -U u +ir 1 /2' 



G5W = ^ + 



e - £ 2 + «r 2 /2 e - E 2 - U 12 + iV 2 /2 ' 

where Uu denotes the Coulomb interaction between two electrons in level 1 and Lq 2 denotes 
that for one electron in level 1 and the other in level 2. To calculate the steady-state electron 
occupation number for the QD, we solve the semiconductor Bloch equations for the two-level 
system coupled to leads and the electromagnetic radiation. [6] We obtain (for Af 2 < #i) 

r\Aq = -ImX - I -[V[f L {e) + rf f R {e)\ImG\{e), 

J 7T 

T 2 iV 2 = ImX- f -\V L 2 f L {e) + T^f R (e)]ImG r 2 (e), 

J 7T 

and 

X(u>) = 2X 2 (N 2 - AM{ l — — , + — -}, 

y ' V n u r -u + iT/2 u r + U 12 -Un -uj + iTjV' 

where u) r = E 2 — E\ is the resonant frequency, r = F 1 + r 2 , and A is the momentum matrix 
element for the inter-level optical transition. 

Note that X{uj) consists of two poles corresponding to optical transitions from initial 
states at E\ and (E 1 + Uu) to final states at E 2 and (E 2 + U\ 2 ), respectively. The two 
transitions have different relative strengths which depend on the averaged steady -state 
occupation number in the ground state, N x . The first term corresponds to the inter- level 
transition of a single electron in the QD [which occurs with a relative probability (1 — Aq)], 
while the second term corresponds to the inter-level transition of a second electron in the QD 
(which occurs with a relative probability Aq) under the influence of the first electron, which 
remains in the ground state at all times. In the latter case, the Coulomb repulsion between 
the two electrons give rise to an energy shift U 12 — Uu. The fact that we have a fractional 



occupancy Aq in a single quantum dot is attributed to the level mixing effect (coupling of 
the QD level to the continuum states in the leads) in the many-body open system. 

The theory is applied to a realistic self-assembled quantum dot (SAQD) device. We 
consider an InAs/InGaAs SAQD system with conical shape. The SAQD is embedded in a 
slab of InGaAs with finite width, W. The slab is then placed in contact with heavily doped 
InGaAs to form an n-i-n structure for infrared detection. Within the effective-mass model [3], 
the QD electron is described by the equation 



m*(p,z) is the position-dependent effective mass, which takes on values of m* G = 0.067m e 
(for GaAs) and mj = 0.024m e (for InAs). The potential V(p, z) is equal to in the InGaAs 
barrier region and Vq inside the InAs QD region. The potential in the depletion layers (which 
separate the slab from the leads) are modelled by an electrostatic potential V d (z) 



For the purpose of constructing the approximate wave functions, we place the system in a 
large cylindrical confining box with length L and radius R (R must be much larger than the 
radius of the cone, r c ). We adopt R = 400A, D = 350 A, V 1 = -0.205eV, and W = 300Afor 
all calculations. We solve the eigen-functions of the effective-mass Hamiltonian via the Ritz 
variational method. The wave functions are expanded in a set of basis functions which are 
chosen to be products of Bessel functions and sine waves 



where k m = mn/L, m=l,2,3... Throughout the paper, the origin of z is set at the middle 
of the confining box. Ji is the Bessel function of order I (I — 0, 1,2, etc.) and a n R is 
the n-th zero of Ji(x). 40 sine functions multiplied by 15 Bessel functions for each angular 
function (I = or 1) are used to diagonalize the Hamiltonian. 



[-v 



1 



V + V(p, z) - eFz]^(p, 0, z) = Eil;(p, 0, z), 



2m*(p, z) 




ipnim(p,<f>,z) = Ji(a n p)e tl4> sm[k m (z + — )] 
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Fig. 1 shows the energy levels of the confined states with / = (solid line) and 1 (dotted 
line) as functions of height h of the QD with base radius fixed at R = 70 A. The other 
material parameters used here are: wetting layer thickness d — 3 A, the conduction-band 
offset Vq = —0.4 eV (this includes the effect of hydrostatic strain due to the lattice mismatch 
between InAs and In0.2Ga0.sAs), and length of the confining box L = 600 A. At least two 
bound states for each angular function (/ = 0, or 1) are found. For infra-red detector 
application, we are seeking an intra-band transition (between the ground and first excited 
state) at an energy around 0.125 eV, which occurs at h = 50 A, for R = 70 A. The tunneling 
rates can be calculated numerically via the stabilization method as described in Ref. 3. 

We only consider the zero temperature and low bias case, where the chemical potential 
at the left leads (//£,) is lower than E 2 , so that the average population in the exited state is 
small. Fig. 2 shows the photocurrent as a function of frequency for various applied voltages: 
solid line (V a = 0.11V), dotted line (V a = 0.12V), and dashed line (V a = 0.13V). The 
parameters used to obtain Fig. 2 are E\ = — 139meV, E 2 = —lAmeV, Uu = lOAmeV, and 
U21 = 7.2meV, which are all calculated based on the effective-mass model. The Fermi level in 
the source and drain region is assumed to E F = 15meV. The broadening of the energy level 
Ei including all tunneling processes (dominated by the acoustic-phonon assisted tunneling 
in this case) is assumed to be Ti = O.OlmeV. The precise value of Ti is not important, 
since photocurrent is not sensitive to Ti. For the excited state, the broadening parameter 
is given by T 2 = T' 2 + T 2 . T 2 is mainly due to radiative and non-radiative recombination 
from interacting with phonons and defects. The actual value depends on the sample quality 
and temperature. Here, we assume T' 2 = lmeV. The other contribution due to the direct 
tunneling is calculated via the stabilization method as described in Ref. [3]. The values 
are found to be = 0A39meV, 0.545meV, and 0.651meV for V a = 0.11,0.12, and 0.13V, 
respectively. The spectrum of photocurrent consists of two peaks centered at frequencies 
ijj = E 2 — Ei, and 00 = E 2 —Ei + Ui 2 — Un. The relative strength of these peaks are determined 
by the average occupation number in the ground state (Ai), which is bias-dependent. As 
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shown in the figure, the Coulomb interaction leads to a double-peak photocurrent spectrum 
with energy separation related to the intra-level and inter-level Coulomb energies (Uu and 

u 12 ). 

For the QDIP characteristics, the photocurrent versus applied bias (J — V curve) is also 
of interest [7-9] . Fig. 3 shows the photocurrent as a function of bias for frequencies at uo r 
(dotted line) and u r + Ui2 — Uu (solid line). Using Eq. (2) for the polarization, we can readily 
understand the behavior of photocurrent. The behavior of the photocurrent is determined 
by the prefactors (1 — Aq)Aq and AqAq for the two poles at uj = u r and uj = uj r + £q 2 — Uu, 
respectively. At very low bias, Aq is small; thus, the magnitude of the solid line is much 
weaker than that of the dotted line. As the applied bias increases, the solid line displays a 
plateau due to the effect of Coulomb blockade on Aq. When the applied bias overcomes the 
charging effect, Aq w 0.5 and the solid line becomes almost identical to the dotted line. 

In this study, we have used tunneling carriers as the source for photocurrent, in contrast 
to the captured carriers typically used in QWIPs. Due to the phonon bottleneck effect [4], 
it is predicted that the capture rate of the electron by the QD will be low. This could 
reduce the performance of QDIPs which use captured carriers as a photocurrent source. 
Using tunneling carriers as the photocurrent source will not have this drawback, and in this 
case we find that the effect of electron correlation leads to a double-peak spectrum and the 
photocurrent is a highly nonlinear function of the carrier density in the ground state Aq. 
Both of these effects must be taken into account in the analysis of photoresponse of QDIPs. 
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Figure Captions 

Fig. 1. Energies of the bound states of a conical InAs/ In , 2 Ga . 8 As QD as functions of 
height h of the QD with base radius fixed at i?o = 70 A. Solid lines: (/ = 0). Dotted lines: 
(1 = 1). 

Fig. 2. Photocurrent as a function of frequency for different applied voltages: V a = 0.1 IV 
(solid line), V a = 0.12V (dotted line) and V a = 0.13V (dashed line). 

Fig. 3. Photocurrent as a function of bias for incident frequencies at uo = E 2 — E\ (dotted 
line) and oo = E 2 — E\ + U\ 2 — Uu (solid line). 
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